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ABSTRACT 


The  fluid  flow  and  forced  convective  heat  transfer  characteristics  of 
a  uniformly  heated  right  circular  cylinder  resting  on  a  flat  plate  in  a 
crossflow  of  air  were  experimentally  investigated. 

The  circumferential  variation  of  the  Nusselt  niomber  on  the  cylinder 
surface  was  obtained  at  Reynolds  numbers  ranging  from  53,000  to  153,000 
using  a  liquid  cirystal  thermographic  technique.   Pressure  coefficients 
were  obtained  at  five-degree  intervals  around  the  cylinder  circumference 
at  Reynolds  numbers  of  105,000  and  156,000.   Flow  visualization  was  con- 
ducted at  a  Reynolds  niimber  of  10,000  in  order  to  qualitatively  examine 
the  flow  field. 

A  qualitative  theory  based  upon  a  compilation  of  previous  results 
reported  for  flow  around  a  "free"  cylinder  and  flow  over  forward  and 
rear-facing  steps  was  developed. 

Experimental  results  obtained  were  in  agreement  with  the  proposed 
theory,  and  heat  transfer  and  pressure  data  correlated  well  with  one 
another.   In  the  region  between  0  and  180  degrees,  the  heat  transfer  and 
pressure  results  correlated  quantitatively  with  the  results  obtained  by 
Meyer  for  a  "free"  cylinder  at  Reynolds  numbers  above  300,000. 
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I.   INTRODUCnON    

The  primary  purpose  of  this  work  was  to  investigate  the  forced  convec- 
tive  heat  transfer  and  fluid  flow  characteristics  of  a  uniformly  heated 
right  circular  cylinder  resting  on  a  flat  plate  in  a  crossflow  of  air. 
This  particular  flow  situation  models  several  situations  of  practical 
interest,  notably,  the  effect  of  wind  on  rocket  motor  storage  containers 
lying  on  the  ground.   In  this  regard,  this  investigation  supports  the 
previous  work  reported  by  Wirzburger  [1]*. 

Local  Nusselt  numbers  were  obtained  at  various  circumferential  loca- 
tions on  the  cylinder  for  Reynolds  numbers  in  the  range  of  53,000  to  153,000 
using  the  liquid  crystal  technique  developed  by  Field  [2] .   A  four-inch 
diameter  cylinder  was  rigidly  attached  to  a  flat  plate  and  placed  in  a  Ic; 
speed  wind  tunnel  with  a  turbulence  intensity  of  approximately  0.5  to  0.7 
percent.   The  test  cylinder  was  constructed  of  a  carbon- impregnated  resis- 
tive paper  having  an  electrical  resistivity  of  approximately  one  ohm-inch. 
The  cylinder  was  coated  with  microencapsulated  cholesteric  liquid  crystals 
and  packed  inside  with  glass  wool  for  the  purposes  of  both  minimizing 
conduction  and  internal  free  convection  losses,  and  providing  resistance 
to  deformation.   The  cylinder  surface  was  uniformly  heated  by  passing  an 
electrical  current  through  the  resistive  paper,  and  data  v;ore  obtainod 
by  observing  the  position  of  the  isotherms  indicated  by  the  liquid  crystals. 

Pressure  coefficients  were  also  obtained  on  the  surface  of  the  cylinder. 
A  rotatable  four-inch  diameter  aluminum  cylinder  was  fitted  with  four  0.025 
inch  static  pressure  taps  90  degrees  apart.   The  cylinder  was  sealed  to 
the  flat  plate  by  meanr.  of  a  Teflon  stri[j  which  <i]  lowed  tht;  cryiindor  to 

♦Numbers  in  brackets  refer  to  references  at  end  of  thesis. 
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rotate,  but  did  not  allow  flow  to  pass  between  the  plate  and  the  cylinder. 
Pressure  data  were  obtained  using  a  firattssassbe^  ^i^nk  referenced  against 
atmospheric  pressure.   Stagnation  pressure  and  test  section  static  pressure 
were  obtained  from  existing  pressure  taps  in  the  wind  tunnel  which  were 
connected  to  the  same  manometer  bank. 

Flow  visualization  was  conducted  in  a  water  tunnel  using  a  1.5-inch 
diameter  plexiglas  cylinder  sealed  against  an  aluminum  flat  plate.   The 
system  was  scaled  for  dynamic  similitude  with  the  heat  transfer  experimental 
apparatus.   The  cylinder  was  equipped  with  three  dye  ports,  two  of  which 
were  located  approximately  30  degrees  either  side  of  the  attachment  point, 
and  the  other  located  at  approximately  the  90  degree  mark.   Due  to  limita- 
tions in  the  water  piimp,  the  maximum  Reynolds  number  which  could  be 
attained  with  this  apparatus  was  approximately  10,000. 

A  qualitative  theory  of  the  flow  field  and  accompanying  heat  transfer 
was  proposed  based  upon  previous  investigations  of  flow  around  a  "free" 
cylinder  combined  with  investigations  of  flow  over  forward  and  rear-facing 
steps. 

The  heat  transfer  and  pressure  data  obtained  were  in  excellent  agree- 
ment with  the  proposed  theory,  and  the  pressure  data  correlated  well  with 
the  heat  transfer  data.   The  flow  visualization  showed  agreement  with 
portions  of  the  proposed  theory,  but  did  not  support  the  heat  transfer 
and  pressure  data  on  certain  regions  of  the  cylinder.   This  is  most 
probably  because  of  the  wide  difference  in  the  Reynolds  number  between 
these  experiments. 


II.   THEORY 


A  thorough  search  of  the  existing  literature  indicated  that  this 
particular  flow  situation  had  not  been  examined.   The  study  of  flow  around 
a  "free"  cylinder  has  been  studied  extensively,  and  an  excellent  bibliog- 
raphy on  the  subject  has  been  presented  by  Meyer  [3] .   Subsequent  to 
Meyer's  work,  Field  [2]  investigated  both  subcritical  and  critical  flow 
around  a  right  circular  cylinder  in  crossflow  using  the  liquid  crystal 
technique  employed  in  this  investigation.   A  brief  description  of  these 
two  flow  situations  is  presented  below. 

In  siabcritical  flow,  a  laminar  boundary  layer  develops  at  the  forward 
stagnation  point  of  the  cylinder  and  grows  to  a  maximum  thickness  at  an 
^ngular  location  of  80  to  85  degrees.   Concurrently,  the  heat  transfer 
coefficient  decreases  as  the  boundary  layer  thickness  increases  in  this 
region.   At  80  to  85  degrees,  the  kinetic  energy  of  the  moving  fluid  is 
no  longer  sufficient  to  overcome  the  adverse  pressure  gradient  present 
on  the  cylinder's  surface,  and  the  laminar  boundary  layer  separates. 
Upon  separation,  a  wake  is  formed  in  the  rear  portion  of  the  cylinder,  and 
the  turbulent  action  in  the  wake  tends  to  circulate  cooler  fluid  to  the 
surface  of  the  cylinder.   As  a  result  of  this  action,  the  heat  transfer 
coefficient  increases  in  this  area. 

In  critical  flow,  which  occurs  at  Reynolds  numbers  above  approximately 
120,000,  the  laminar  boundary  layer  again  grows  from  the  forward  stagnation 
point  to  the  80  to  85  degree  point  whore  separation  again  occurs.   However, 
in  this  case  the  kinetic  energy  of  the  fluid  is  enough  to  overcome  the 
adverse  pressure  gradient,  and  the  flow  reattaches  itself  to  the  cyliiuujr 


10 


surface  at  a  point  approximately  15  degrees  aft  of  the  laminar  separation 
point.   This  separation  "bii)ble~'  Is  'tii52i3-46!£:r-i3i^^        highly  turbulent 
mixing  and  results  in  a  sudden  increase  in  heat  transfer.   From  the  point 
of  reattachment,  a  turbulent  boundary  layer  begins  to  grow  to  a  point  in 
the  region  of  110  to  130  degrees  where  it  separates.   The  heat  transfer 
coefficient  decreases  as  the  turbulent  boundary  layer  thickens.   Upon 
separation  of  the  turbulent  boundary  layer,  the  action  of  the  wake  causes 
an  increase  in  the  heat  transfer  coefficient. 

It  should  be  noted  that  the  above  description  is  confined  to  a  right 
circular  cylinder  which  is  freely  suspended  in  crossflow  such  that  flow 
can  pass  over  both  sides.   Hence,  the  flow  pattern  is  spacially  symmetrical. 
There  was  no  reason  to  expect  this  to  be  the  case  in  this  investigation, 
however,  and  an  attempt  was  made  to  gain  some  understanding  of  what  results 
might  be  expected  in  the  regions  of  the  cylinder  near  the  flat  plate.   This 
was  accomplished  by  examining  the  cases  of  flow  over  forward- facing  and 
rear-facing  steps. 

The  problem  of  flow  over  forward-facing  and  rear-facing  steps  has  been 
examined  in  various  lights.   Several  investigators  have  focused  on  the 
determination  of  the  location  of  the  separation  points  and  reattachment 
points  [4-9] .   Others  have  examined  the  heat  transfer  and  pressure  dis- 
tribution v;ithin  the  separated  regions  [10-17] .   However,  these  investiga- 
tions were   primarily  concerned  with  the  regions  of  the  flat  plate  either 
forward  or  aft  of  the  step  and  did  not  deal  with  the  step  itself.   Since  a 
cylinder  resting  on  flat  surface  represents  a  type  of  "step",  this  region 
was  of  most  interest  in  this  work.   A  few  investigators  did  conduct  flow 
visualization  experiments  in  conjunction  wi,th  tlioir  other  goals,  however, 
and  their  results  proved  to  be  of  benefit  in  predicting  the  naLure  of  the 
flow  in  the  regions  of  the  cylinder  near  the  flat  plate. 
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Of  particular  interest  were  the  findings  of  Abbott  and  Kline  [5]  for 
rear- facing  steps,  and  those  of  Luznaii^teai,'  '^asS'^SsoScfesev  [10]  for  forward- 
facing  steps.   The  latter  investigated  the  zone  of  turbulent  boundary  layer 
separation  ahead  of  a  forward- facing  step  at  Reynolds  numbers  of  17,000  to 
420,000  based  upon  step  height.   Although  their  visualization  technique  was 
not  described,  their  experiments  were  conducted  in  air.   The  results  of 
this  investigation  are  shown  in  Figure  1.   It  should  be  noted  that,  in 
addition  to  the  normal  separation  vortex  (points  1-2-1) ,  another  small 
"trapped"  vortex  was  noted  in  the  region  between  points  3  and  4. 

Abbott  and  Kline  conducted  flow  visualization  experiments  using  the  dye 
injection  method  on  a  variable  height  backward- facing  step  placed  in  a 
water  tunnel.   The  boundary  layer  at  separation  from  the  step  was  turbulent 
as  is  the  case  for  critical  flow  around  a  cylinder.   The  experiments  were 
conducted  in  a  Reynolds  number  range  of  20,000  to  50,000  based  upon  the 
height  of  the  water  channel  above  the  step.   Their  results  indicated  that 
three  regions  of  flow  existed  downstream  of  the  step.   Region  I  was 
characterized  by  one  or  more  vortices  rotating  about  an  axis  normal  to  the 
vertical  wall  and  parallel  to  the  floor.   These  vortices  were  three- 
dimensional  in  nature,  and  adjacent  vortices  were  counter-rotating  and  not 
necessarily  the  same  size.   Region  II  contained  the  classical  pattern  of 
separation,  i.e.  flow  moving  two-dimensionally  upstream  next  to  the  wall 
and  downstrcani  adjacent  to  the  through  flow.   Region  III  was  of  no  interest 
to  the  present  investigation  because  it  did  not  directly  affect  the  step. 
The  three  regions  are  depicted  both  in  plan  and  elevation  views  in  Figure 
2. 

From  a  compilation  of  the  results  described  above  and  those  observed 
for  a  "free"  cylinder,  the  following  theory  was  developed.  ' 
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Figure  1.   Flow  pattern  ahead  of  a  step  as  observed  by 
Luzhanskiy  and  Solntsev. 
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Figuro  2.   Flow  pattern  at  the  rear  of  a  step  as  observed 
by  /Ujbot'.t  and  Kl  inc. 
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As  the  flow  proceeds  along  the  flat  plate,  a  boundary  layer  will  develop 
and  grow  on  the  plate  until,  at  some  point  upstream  of  the  cylinder,  this 
boundary  layer  will  separate.   At  the  flow  velocities  to  be  used  in  the 
experimental  portion  of  this  investigation,  this  boundary  layer  can  be 
anticipated  to  be  laminar.   From  this  point,  it  is  theorized  that  the 
flow  will  be  characterized  by  a  behavior  similar  to  that  observed  by 
Luzhanskiy  and  Solntsev,  i.e.  a  portion  of  the  flow  will  reattach  to  the 
cylinder  at  some  point  on  the  forward  portion  of  its  circumference,  and 
the  remainder  of  the  flow  will  assume  the  classical  separation  pattern  with 
the  addition  of  a  system  of  "trapped"  vorticies  near  the  intersection  of 
the  cylinder  and  the  flat  plate.   The  flow  which  reattaches  itself  to  the 
cylinder  then  will  assume  the  classical  pattern  described  previously  for 
flow  around  a  "free"  cylinder,  with  the  differentiation  between  subcritical 
and  critical  flow  being  dependent  upon  the  Reynolds  number.   Regardless  of 
whether  the  flow  is  subcritical  or  critical,  it  will  separate.   Upon 
separation,  it  is  assumed  that  a  pattern  similar  to  that  observed  by 
Abbott  and  Kline  will  be  present. 

From  the  above  theory  of  the  flow  phenomenon,  the  heat  transfer  results 
can  be  predicted.   For  ease  of  explanation  the  cylinder  will  be  described 
as  consisting  of  three  regions.   Region  I  consists  of  that  portion  of  the 
cylinder  circuinforence  betv;cen  the  point  of  attachment  of  the  cylinder  to 
the  flat  plate  and  the  point  at  which  the  separated  flow  from  the  plate 
reattaches  to  the  cylinder.   This  latter  point  will  hereafter  be  referred 
to  as  the  stagnation  point  for  simplicity.   Region  II  consists  of  that 
portion  of  the  cylinder  circuinlerence  boLwoen  the  stagnation  point  and  the 
point  of  final  boundary  layer  separation.   Region  III  consists  of  that 
portion  of  the  cylinder  circumference  aft  of  the  final  separation  point. 
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The  phenomenon  which  is  proposed  to  occur  in  Region  I  is  depicted  in 
Figure  3.   Since  all  the  fluid  \rj^xljLii'^^ii^its:*<i;&s^,tMjJs:>n   zone  is  "trapped," 
heat  liberated  from  the  cylinder  surface  must  be  conducted  from  vortex  to 
vortex  within  the  separation  region  and  ultimately  conducted  to  the  free 
stream  through  the  free  shear  layer  separating  the  separation  zone  from 
the  free  stream.   This  being  the  case,  the  bulk  temperature  in  the  larger 
vortices  must  be  less  than  the  bulk  temperature  in  the  smaller  vortices 
in  order  for  a  driving  potential  for  conduction  to  exist.   Convection, 
however,  also  plays  a  role  in  this  process  in  that  a  redistribution  of 
energy  within  individual  vortices  is  accomplished  by  convection.   It  is 
assumed  that  the  vortices  will  become  smaller  and  weaker  as  the  attachment 
point  is  approached  such  that,  very  near  this  point,  the  heat  transfer  will 
approach  a  pure  conduction  case.   Hence,  the  resistance  to  heat  transfer 
should  increase  from  the  stagnation  point  to  the  attachment  point.   If  the 
thermal  resistance  increases,  the  wall  temperature  must  also  increase  in 
order  to  insure  that  a  constant  heat  flux  is  transferred  to  the  fluid.   If 
the  wall  temperature  increases  and  the  free  stream  temperature  and  surface 
heat  flux  remain  constant,  the  heat  transfer  coefficient  must  decrease. 
Hence,  in  the  absence  of  other  effects,  the  heat  transfer  coefficient 
should  decrease  from  a  relative  maximum  at  the  stagnation  point  to  a  zero 
value  at  the  attachment  point  (assuming  the  flat  plate  is  adiabatic) . 
Superimposed  upon  on  this  phenomenon,  hov/ever,  is  the  behavior  of  the  heat 
transfer  coefficient  near  the  intersection  of  two  vortices.   The  region 
between  points  1  and  3  of  Figure  3  is  illustrative  of  this  behavior.   At 
point  1,  the  stagnation  point,  the  fluid  in  the  vortex  which  impacts  on 
the  cylinder  surface  is  relatively  cool  because  some  of  its  energy  has 
just  been  conducted  to  the  free  stream.   As  the  fluid  proceeds  toward 
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point  2,  it  is  heated  and  its  thermal  resistance  is  increased;  hence,  the 
heat  transfer  coefficient  is  decreasea' jDsV'^^tiaas'SEgion.   Since  adjacent 
vortices  are  counter-rotating,  a  similar  behavior  should  occur  between 
points  3  and  2.   Therefore,  point  2  should  appear  as  a  relative  minimum 
in  the  heat  transfer  coefficient,  and  the  heat  transfer  coefficient  should 
then  increase  to  a  relative  maximiam  at  point  3.   At  point  3,  the  cycle 
described  above  should  begin  again,  this  time  involving  points  3,  4,  and  5. 
It  should  be  noted  that,  because  the  vortices  become  smaller  and  weaker 
as  the  point  of  attachment  is  approached,  this  secondary  effect  may  become 
experimentally  imperceptible  in  the  region  near  this  point.   It  should  also 
be  noted  that  the  number  of  vortices  depicted  in  the  schematic  does  not 
consititute  a  prediction  of  how  many  will  actually  be  present. 

Region  II  is  characterized  by  the  results  reported  by  Field  [2]  for 
subcritical  or  critical  flow,  whichever  may  be  the  case. 

Region  III  is  theorized  to  indicate  an  initial  increase  in  the  heat 
transfer  coefficient  due  to  the  turbulent  action  of  the  wake.   However, 
at  some  point,  this  increase  will  be  overridden  by  the  phenomenon  described 
for  Region  I,  and  the  heat  transfer  coefficient  will  decrease  to  a  zero 
value  at  the  attachment  Doint. 

Schematics  depicting  the  theorized  flow  situation  and  accompanying  heat 
transfer  data  for  both  subcritical  and  critical  flows  are  s'.hown  as  Figures 
4  and  5  respectively.  Tnc    locatioiis  at  which  various  plienomenon  occur 
are  conjectural  and  should  not  be  construed  as  being  exact  predicted 
values. 
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rigiiro  5.   Proposed  flow  pattern  for  critical  RcynoldG  numborf; 
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III.   EXPERIMENTAL  APPARATUS 

A.   HEAT  TRANSFER  EXPERIMENT 

1.  Liquid  Crystals 

Microencapsulated  cholesteric  liquid  crystals  were  used  in  this 
investigation.   This  substance,  when  applied  to  a  dark  background,  provides 
a  visual  display  of  temperature  by  means  of  a  change  of  color.   The  surface 

coated  with  liquid  crystals  traverses  the  visible  color  spectrum  within  a 

o 
certain  temperature  band  with  an  accuracy  of  approximately  0.2  F.   The 

size  of  this  band  is  deteirmined  by  the  formulation  of  the  cholesterol  esters 

in  the  mixture.   Only  one  liquid  crystal  formulation  was  used  in  this 

investigation,  NCR  designation  S-43.   The  batch  used  had  been  calibrated 

previously  by  Field  [2] ,  and  the  results  of  this  calibration  indicated 

o 
that  liquid  crystal  S-43  transitioned  to  a  red  color  at  109.8  F,  to  a 

o  o 

green  color  at  111.0  F,  and  to  a  blue  color  at  112.3  F.   An  excellent  list 

of  references  on  the  theory  and  use  of  liquid  crystals  can  be  found  in 

Field's  work. 

2.  Cylinder  Design  and  Construction 

The  cylinder  utilized  in  the  collection  of  heat  transfer  data  was 
similar  in  design  to  that  used  by  Field.   The  actual  working  section  was 
constructed  from  a  12-inch  square  sheet  of  Armstrong  Temsheet,  a  carbon- 
impregnated,  electrically  resistive  paper  having  an  electrical  resistance 
of  approximately  25  ohms  per  square  and  a  nominal  thickness  of  0.039  inches. 
The  heat  generated  when  a  constant  current  is  passed  through  the  paper  is 
uniform  to  within  two  percent  from  point  to  point  over  a  large  area.   This 
sheet  was  fasluoned  into  a  cylinder  having  an  outside  diameter  of  4.03 
inches  and  a  height  of  12  inches. 
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Aluminum  conducting  tape  was  applied  as  electrodes  circumferential ly 
on  the  inner  surface  of  this  cylinder  to  insure  uniform  application  of  volt- 
age.  These  electrodes  were  located  three  inches  from  the  top  and  bottom  of 
the  cylinder  in  order  that  edge  conduction  effects  would  be  eliminated. 
This  established  a  six  inch  height  of  cylinder  in  which  data  would  be  taken. 
Silver  conducting  paint  was  applied  to  the  interface  between  the  tape  and 
the  paper  to  insure  intimate  electrical  contact  was  maintained  between  the 
two  surfaces. 

The  base  assembly  consisted  of  two  4.03-inch  diameter  wooden 
cylinders,  each  14.75  inches  long,  rigidly  attached  to  each  other  by  an 
eight-inch  long  wooden  bar.   The  outside  surface  of  this  bar  was  of  the 
same  curvature  as  the  cylinders  in  order  to  insure  that  no  discontinuities 
existed  near  the  point  of  attachment  to  the  flat  plate.   The  purpose  of 
this  bar  was  to  maintain  alignment  between  the  two  base  cylinders,  and  to 
allow  the  Temsheet  cylinder  to  be  taped  along  its  length  to  the  bar .   The 
entire  assembly  was  slotted  to  accept  a  1/4-inch  square  plexiglas  rod  along 
its  length.   This  rod  accepted  screws  from  the  flat  plate  to  insure  that 
an  acceptable  seal  was  maintained  at  the  point  of  attachm.ent.   The  base 
assembly  is  shoum  in  Figure  6. 

Each  base  cylinder  v/as  of  sufficient  length  such  that  it  would 
extend  through  the  floor/ceiling  of  tho  v/ind  tunnel  where  it  could  be 
clamped  thus  providing  added  rigidity  to  the  assembly.   The  other  ends  of 
the  base  cylinders  were  turned  down  approximately  0.04  inches  over  a  two- 
inch  length.   This  allowed  the  paper  to  be  attached  to  the  bases  with 
double-backed  tape  and  present  a  constant  diameter  of  4.03  inches  along 
the  entire  length  of  tho  assembly.   The  bottom  base  was  drilled  througli 
its  axis  v/ith  a  one-inch  diameter  hole  through  which  the  electrical  leads 
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Figure  6.   Temsheet  cylinder  base  assembly 
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could  be  passed  to  the  power  supply;  and  the  top  base  was  drilled  through 
its  axis  with  a  two-inch  diameter  hole  to  provide  a  means  of  filling  the 
> Temsheet  cylinder  with  glass  wool.   The  glass  wool  served  the  dual  purpose 
of  providing  insulation  against  possible  internal  heat  losses  due  to  con- 
duction or  natural  convection  and  of  providing  the  paper  cylinder  with 
added  resistance  to  deformation. 

Two  coats  of  slurry-based  liquid  crystal  S-43  were  applied  to  the 
Temsheet  cylinder,  and  it  was  marked  at  five-degree  intervals  for  accurate 
location  of  the  liquid  cirystal  isotherms  during  the  experiment. 

The  flat  plate  was  constructed  from  a  sheet  of  1/2-inch  plexiglas 
and  was  36  inches  long  and  32  inches  high.   Plexiglas  was  selected  for  this 
constructipn  in  order  that  the  entire  circumference  of  the  cylinder  could 
be  viewed  from  outside  the  wind  tunnel.   The  leading  edge  of  the  plate  was 
beveled  approximately  14  degrees  to  insure  that  a  sharp  leading  edge  v.'as 
presented  to  the  flow.   The  point  of  attachment  of  the  cylinder  was  located 
on  a  vertical  line  18  inches  downstream  of  the  leading  edge.   At  this 
point,  the  plate  was  notched  along  its  height  to  accept  the  plexiglas  rod 
on  the  cylinder  assembly,  and  attachment  was  accomplished  by  four  screws. 
This  design  insured  that  flow  could  not  pass  between  the  cylinder  and  the 
,   plate. 

The  wind  tunnel  in  which  this  cylinder  assembly  was  placed  was  the 
low  speed  Aerolab  wind  tunnel  located  in  the  Aeronautics  Laboratory  in 
Halligan  Hall,  U.  S.  Naval  Postgraduate  School.   This  tunnel  is  powered  by 
a  100  hp  electric  motor  and  is  capable  of  producing  wind  velocities  of  up 
to  200  mph  v/ith  a  clear  test  section.   The  tunnel  has  a  turbulence  level 
as  measured  by  hot  wire  anemomotcr  of  api)roximatcly  0.5  to  0.7  percent.   A 
schematic  of  the  wind  tunnel  v;ith  pertinent  dimensions  is  shown  in  Figure  7. 
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The  combined  cylinder  and  plate  assembly  was  located  in  the  wind 
tunnel  test  section  such  that  the  longitudinal  axis  of  the  flat  plate  lay 
along  a  line  eight  inches  from  one  wall  of  the  test  section.   By  off- 
'  centering  the  assembly  in  this  manner,  the  effects  of  the  flow  along  the 
wind  tunnel  wall  on  the  cylinder  were  effectively  eliminated. 
B.   PRESSURE  EXPERIMENT 

In  order  to  obtain  static  pressure  measurements  of  this  particular 
flow  situation,  a  4.03-inch  diameter  cylinder  was  constructed  of  aluminum. 
This  cylinder  was  fitted  with  four  0.025-inch  surface  pressure  taps  spaced 
at  90-degree  intervals  around  the  circumference  of  the  cylinder  at  approxi- 
mately mid-height  of  the  wind  tunnel  test  section.   In  order  that  this 
cylinder  could  be  rotated  and  a  seal  still  be  maintained,  a  Teflon  strip 
was  inserted  into  the  notch  in  the  flat  plate  described  previously.   This 
strip  provided  for  a  snug  fit  between  the  cylinder  and  the  plate,  and  the 
flexibility  of  the  Teflon  allowed  for  relatively  easy  rotation.   The 
cylinder  was  scribed  as  five-degree  inter-vals  about  its  circumference  at 
the  point  where  it  intersected  the  wind  tunnel  test  section  floor. 
Reference  lines  were  scribed  on  the  floor  at  90-degree  intervals  so  that 
the  matching  of  the  scribe  marks  on  the  cylinder  with  these  lines  provided 
a  means  of  accurately  establishing  the  angular  location  of  the  pressure 
taps.   The  prer.surc  cylinder  is  shown  in  Figure  8. 

The  surface  pressure  taps  were  connected  to  a  water  manometer  bank  by 
means  of  plastic  tubing.   The  monometer  bank  was  referenced  against 
atmospheric  pressure.   Wind  tunnel  plenum  chamber  static  pressure  and  test 
section  static  pressure  were  measured  on  the  same  manometer  bank  through 
pre-existing  connections  on  the  wind  tunnel. 
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Figure  8.   Cylinder  used  to  collect  pressure  data 
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C.   FLOW  VISUALIZATION  EXPERIMENT 

The  apparatus  used  for  flow  visualization  consisted  of  a  converging 
wooden  water  channel  15  feet  in  length  and  13  inches  deep  with  a  plexiglas 
test  section  13  inches  long,  12  inches  wide,  and  13  inches  deep.   A  hollow, 
1.5-inch  diameter  plexiglas  cylinder  was  attached  between  the  test  section 
walls  with  its  longitudinal  axis  located  eight  inches  from  the  test  section 
floor.   In  order  to  preserve  dynamic  similitude,  a  1/4-inch  thick  steel 
flat  plate  was  constructed  to  the  scale  of  the  ratio  of  the  water  tunnel 
cylinder  diameter  to  the  wind  tunnel  cylinder  diameter  (0.375:1).   Thus, 
this  plate  was  13.5  inches  long  and  12  inches  wide.   The  leading  edge  was 
beveled  to  a  14  degree  angle  as  was  the  larger  flat  plate  used  in  the  wind 
tunnel.   The  flat  plate  was  placed  atop  the  cylinder,  and  the  point  of 
attachment  was  sealed  by  a  layer  of  tape.   The  cylinder  was  fitted  with 
three  dye  ports  which  were  connected  by  rubber  tubing  to  three  separate 
dye  reservoirs  with  individual  valving  located  approximately  three  feet 
higher  than  the  cylinder.   This  additional  height  insured  that  a  constant 
head  was  maintained  on  the  injection  system,  and  the  individual  valving  was 
necessary  to  allow  each  port  to  be  operated  separately  in  order  that  the 
flow  field  at  various  sections  of  the  cylinder  could  be  observed  independently. 

VVatcr  was  pumped  into  the  channel  by  means  of  an  irrigation  pump  directly 
coupled  to  a  VJestinghouse  Type  SK  dynamometer  running  as  a  motor.   With  the 
motor  running  at  1500  RPM,  the  pump  output  was  approximately  1000  gallons 
per  minute.   The  depth  of  water  in  the  channel  was  controlled  by  a  tilting 
weir  at  the  end  of  the  channel.   For  this  particular  experiment,  the  water 
depth  was  maintained  at  10  inches.   This  allowed  the  entire  apparatus  to 
be  immersed  to  a  minimum  deptli  in  order  that  free  surface  effects  could  be 
minimized. 
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IV.   EXPERIMENTAL  PROCEDURES 

A.   HEAT  TRANSFER  EXPERIMENT 

Prior  to  conducting  the  experiment,  the  thermocouple  used  to  measure 
airstream  temperature  and  its  ice  bath  reference  junction  were  examined, 
and  the  micromanometer  used  to  measure  wind  velocity  was  leveled  and  care- 
fully zeroed.   The  test  cylinder  power  leads  were  then  connected  to  a 
Weston  ohmmeter,  and  the  cylinder  resistance  was  determined.   The  leads 
were  then  switched  to  a  Lambda  regulated  power  supply  model  LK345A,  and 
voltage  was  applied  until  the  temperature  of  the  entire  cylinder  surface 
fell  within  the  event  temperature  range  of  liquid  crystal  S-43  causing 
the  cylinder  to  undergo  a  color  change.   This  procedure  was  necessary  to 
insure  that  the  glass  v.'ool  packing  v.'as  adjusted  properly  so  that  regions 
of  temperature  discontinuity  did  not  exist.   The  wind  tunnel  was  then 
started  and  brought  to  a  predetermined  velocity  at  which  it  was  maintained 
throughout  the  run.   Applied  voltage  was  adjusted  until  distinct  positions 
of  temperature  transition  could  be  observed.   After  steady  state  was 
reached,  the  angular  locations  and  magnitudes  of  these  isotherms  were 
recorded  in  addition  to  che  values  of  applied  voltage  and  airstream  tempera- 
ture.  A  steady  state  condition  was  usually  attained  in  from  five  to  ten 
minutes  dependent  uron  the  wind  velocity.   The  voltage  was  then  changed 
until  the  angular  location  of  the  isotherms  had  changed  a  significant 
amount,  and  data  were  recorded  again.   This  process  was  repeated  until 
isotherms  had  been  observed  at  reasonable  intervals  around  the  circumference 
of  the  cylinder.   Through  this  procedure,  a  plot  of  local  Nusselt  number 
versus  angular  location  could  be  quickly  and  easily  obtained. 
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Runs  were  made  at  nine  Reynolds  numbers  ranging  between  53,000  and 
153,000.   Each  riin  was  accomplished  in  from  60  to  90  minutes,  this  time 
being  dependent  upon  Reynolds  number.   For  all  runs  the  points  of  tempera- 
ture transition  were  easily  distinguishable  around  the  entire  circumference 
of  the  cylinder  except  in  the  region  from  0  to  45  degrees  and  in  the  region 
from  180  to  210  degrees.   In  these  regions,  the  temperature  gradient  was 
very  shallow,  and  the  entire  region  appeared  isothermal.   In  addition, 
data  were  not  obtainable  in  the  regions  15  degrees  either  side  of  the 
attachment  point  (the  270  degree  mark)  because  a  break  in  the  resistive 
paper  was  necessary  in  these  regions  to  facilitate  the  attachment  of  the 
cylinder  to  the  flat  plate. 

B.   PRESSURE  EXPERIMENT 

Pressure  data  v;ere  obtained  at  Reynolds  numbers  of  105,000  and  156,000. 
Initially,  the  manometer  bank  was  leveled,  all  leads  were  checked  for 
tightness,  and  the  air  temperature  was  recorded  for  evaluation  of  transport 
properties.   The  pressure  taps  were  initially  located  at  0,  90,  180,  and 
270  degrees.   The  wind  tunnel  was  started,  and  the  velocity  was  increased 
to  a  predetermined  value.   Values  were  recorded  for  each  pressure  tap  on 
the  cylinder,  the  wind  tunnel  pieiiuia  chairiber  static  pressure,  the  test 
section  static  pressure,  and  atmospheric  pressure.   The  cylinder  was  then 
rotated  in  five-degree  increments,  with  data  being  recorded  at  each 
increment,  until  the  entire  circumference  of  the  cylinder  had  been  traversed. 
Data  were  then  reduced,  and  a  preliminary  plot  of  pressure  coefficient 
versus  angular  location  v.'as  made.   This  plot  v/as  examined,  and  regions 
v/arrantiny  closer  investigation  v/ere  identified.   Data  were  then  f>btained 
at  2  1/2-dcgrcc  increments  within  tlieso  regions. 
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C.   FLOW  VISUALIZATION  EXPERIMENT 

The  motor  was  started  and  brought  up  to  a  speed  of  1500  RPM.   The 
tilting  weir  was  adjusted  such  that  the  height  of  the  water  in  the  channel 
was  10  inches,  and  the  flow  was  allowed  to  steady  out.   The  dye  was  then 
injected  through  the  existing  ports  in  a  trial  and  error  fashion  until  the 
flow  could  be  visualized  clearly.   Most  regions  of  interest  could  be 
observed  in  this  manner.   In  those  regions  which  required  closer  observa- 
tion, dye  was  injected  directly  into  the  flow  by  means  of  a  hypodermic 
needle. 

It  became  necessary  to  interrupt  the  procedure  periodically  due  to  the 
presence  of  large  air  bubbles  which  coalesced  on  the  flat  plate  downstream 
of  the  cylinder.   These  bubbles  were  thought  to  have  been  caused  by 
entrained  air  in  the  water  which  left  solution  in  regions  of  low  press\ire. 
These  bubbles  were  swept  away  by  hand  before  the  experiment  was  continued. 
It  should  be  noted  that  the  effect  of  these  air  bubbles  was  not  completely 
deleterious.   In  certain  regions  on  the  cylinder  surface  where  the  dye  did 
not  provide  adequate  information  regarding  the  existing  flow  conditions, 
the  bxobbles  could  be  used  as  a  flow  visualization  aid.   The  location  of 
the  laminar  separation  point  was  found  in  this  manner. 


30 


V.   RESULTS 

A.   HEAT  TRANSFER  EXPERIMENT 

Raw  data  and  sample  calculations  from  the  heat  transfer  experiments 
are  contained  in  APPendix  A.   This  section  contains  graphical  results  and 
comparisons  of  the  local  Nusselt  number  versus  angular  location  for  various 
Reynolds  numbers,  and  of  the  average  Nusselt  number  versus  Reynolds  number. 

Figure  9  shows  the  results  of  the  heat  transfer  experiment  for  Reynolds 
numbers  of  53,000,  103,500,  and  153,000.   The  ratio  of  boundary  layer 
thickness  on  the  flat  plate  to  cylinder  diameter  was  approximately  0.3. 
Since  the  shape  of  these  curves  did  not  vary  throughout  the  Reynolds  num- 
ber range  tested,  the  curves  presented  are  considered  representative. 
There  are  several  points  of  interest  to  be  noted  from  these  curves. 

First,  the  flow  around  the  cylinder  exhibited  critical  behavior  at  all 
Reynolds  numbers  tested.   Otherwise,  there  was  good  agreement  between 
these  results  and  the  proposed  theory.   The  stagnation  point  was  located 
in  the  region  between  5  and  20  degrees  and  appeared  to  move  toward  0  degrees 
with  increasing  Reynolds  number.   The  point  of  laminar  separation  was 
relatively  stationary  at  approximately  75  degrees.   Reattachment  occurred 
consistently  at  approximately  92  degrees,  and  final  separation  occurred 
consistently  at  approximately  140  degrees.   As  Reynolds  numbor  increased, 
the  point  of  maximum  heat  transfer  shifted  from  the  stagnation  point  to 
the  point  of  reattachment  of  the  turbulent  boundary  layer.   Of  particular 
note  are  the  tv/o  points  of  relative  maxima  ]ocated  at  approximately  240 
and  300  degrees  respectively.   Since  the  point  of  attachment  of  the 
cylinder  to  the  flat  plate  was  Laken  to  bo  the  270  degree  mark,  these 
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points  lie  in  the  regions  in  which  this  type  of  behavior  was  postulated 
by  theory. 

Figures  10  and  11  are  a  comparison  of  the  heat  transfer  data  obtained 
in  this  investigation  with  that  obtained  by  Field  [2]  at  two  nearly  equal 
Reynolds  numbers ,  one  in  the  subcritical  range  and  one  in  the  critical 
range.   Figure  11  also  contains  a  plot  of  the  data  obtained  by  Meyer  [3] 
at  a  Reynolds  number  of  495,000.   Additionally,  a  comparison  of  average 
Nusselt  number  versus  Reynolds  number  for  this  investigation  and  Field's 
investigation  is  shown  as  Figure  12.   The  average  Nusselt  number  was 
computed  from  the  formula 


_1    /"2Tr 

'-    jo 


Nu  =  —   /     Nu  de 


which  was  numerically  integrated  by  computer  using  Simpson's  Rule  of 
Quadrature  with  a  step  size  of  5  degrees. 

It  is  readily  appail'ent  that  the  presence  of  the  flat  plate  radically 
altered  the  local  heat  transfer  coefficient  around  the  cylinder  especially 
in  the  classical  subcritical  Reynolds  number  range  (Reynolds  number  less 
than  120,000).   In  the  critical  Reynolds  number  range,  the  overall  trend 
in  the  curves  was  similar,  but  the  comparative  values  of  the  Nusselt 
number  at  locations  of  interest  were  widely  differing.   In  the  region 
from  0  to  120  degrees,  the  values  found  in  this  investigation  for  Reynolds 
numbers  in  the  range  between  53,000  and  153,000  were  on  the  order  of  those 
found  by  Meyer  for  a  "free"  cylinder  at  Reynolds  numbers  above  300,000. 
Field's  point  of  maximum  heat  transfer  occurred  at  the  180  degree  mark 
whereas,  this  investigation  showed  the  maximum  value  to  occur  at  the 
point  of  reattachment  of  the  turbulent  boundary  layer.   Meyer  also  found 
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the  maximum  value  to  occur  at  the  reattachment  point,  but  this  phenomenon 
was  also  observed  at  Reynolds  numbers  at)ove"j(t)b,u00. 

It  is  also  interesting  to  note  that,  despite  the  lack  of  correlation 
between  this  data  and  Field's  data,  the  plots  of  average  Nusselt  number 
versus  Reynolds  number  have  nearly  the  same  slope  in  the  Reynolds  number 
range  between  50,000  and  150,000,  with  Field's  average  values  being 
higher  than  those  obtained  in  this  investigation  at  all  Reynolds  numbers. 

B.   PRESSURE  EXPERIMENT 

Raw  data  and  sample  calculations  from  the  two  pressure  runs  are  con- 
tained in  Appendix  B.   This  section  contains  graphical  results  of  the 
pressure  coefficient  versus  angular  location  for  Reynolds  numbers  of 
105,000  and  156,000  (Figure  13) ,  and  a  comparison  of  the  pressure 
coefficient  curve  at  a  Reynolds  number  of  155,000  with  rhose  obtained 
by  Meyer  [3]  at  Reynolds  numbers  of  153,000  and  495,000  (Figure  14). 

The  pressure  distributions  at  the  two  Reynolds  numbers  tested  were  in 
good  agreement  with  the  heat  transfer  results  obtained  at  similar  Reynolds 
numbers.   This  is  especially  evident  in  the  higher  Reynolds  niamber  case. 
For  this  case,  the  pressure  coefficient  reached  a  minimum  at  73  degrees, 
and  the  subsequent  upward  slope  contained  a  double  inflection  character- 
istic of  the  separation  "bubble"  region  between  75  and  85  degrees.   The 
first  inflection  point  occurred  at  approximately  78  degrees  which  corre- 
sponds precisely  to  the  laminar  separation  point  found  in  the  heat  transfer 
experiment.   Likewise,  the  second  inflection  point  occurred  at  approximately 
90  degrees  which  corresponds  closely  with  the  92-degrGe  location  of  the 
reattachment  point  found  in  the  heat  transfer  experiment.   Additionally, 
the  pressure  coefficient  curve  began  to  flatten  at  a£ij)roximately  135 
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degrees  which  corresponds  favorably  with  the  point  of  final  separation 
observed  at  140  degrees  in  the  heat  t.ix.iii^42x- 5£&;t>eriinent .   Finally,  fluctua- 
tions in  the  pressure  coefficient  curves  were  noted  in  the  regions  cor- 
responding to  the  locations  of  the  "trapped"  vortices  which  reinforces 
the  heat  transfer  data  in  this  region. 

A  comparison  of  the  pressure  coefficient  curves  obtained  at  a  Reynolds 
nijinber  of  156,000  with  those  of  Meyer  obtained  on  one  side  of  a  "free" 
cylinder  at  Reynolds  numbers  of  153,000  and  495,000  indicated  that  the 
behavior  exhibited  by  the  cylinder  attached  to  the  flat  plate  more  closely 
patterned  that  of  a  "free"  cylinder  at  much  higher  Reynolds  numbers.   In 
the  region  between  0  and  180  degrees,  the  magnitude  of  the  pressure 
coefficients  obtained  at  a  Reynolds  number  of  156,000  were  nearly  identical 
to  those  of  Meyer  obtained  at  a  Reynolds  number  of  495,000,  but  offset 
approximately  10  degrees  throughout  the  entire  region.   This  correlates 
with  the  fact  that  the  magnitudes  of  the  heat  transfer  coefficient  at  this 
Reynolds  number  followed  closely  those  of  Meyer  at  Reynolds  numbers  above 
300,000. 

C.   FLOW  VISUALIZATION  EXPERIMENT 

Because  of  the  limitation  in  the  pump  capacity,  the  highest  Reynolds 
number  which  could  be  obtained  in  this  experiment  was  approximately  10,000. 
This  was  v/ell  below  the  Reynolds  number  range  examined  in  both  the  heat 
transfer  and  pressure  pxperiments,  and,  as  a  result,  correlation  was 
difficult.   However,  the  behavior  observed  agreed  well  with  the  proposed 
theory  for  the  subcritical  range. 

Referring  to  the  regions  of  the  cylinder  described  in  the  proposed 
theory.  Region  I  behaved  almost  exactly  as  predicted.   The  flow  was 
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observed  to  separate  from  the  flat  plate  and  reattach  to  the  cylinder  at 
a  point  just  aft  of  the  0  degree  marl. . -»S*^-H'?:;A^.his  separated  region,  a 
large  turbulent  separation  vortex  was  observed,  and  a  series  of  small, 
circular,  three-dimensional  vortices  were  clearly  evident  at  approximately 
30  to  60  degrees  from  the  point  of  attachment  of  the  cylinder  to  the 
plate. 

Region  II  exhibited  behavior  characteristic  of  subcritical  flow  around 
a  "free"  cylinder.   No  separation  "bubble"  was  noted,  probably  due  to  the 
fact  that  the  Reynolds  number  was  too  low.   However,  laminar  boundary  layer 
separation  was  clearly  observed  at  approximately  80  degrees.   This  behavior 
was  as  predicted  by  the  proposed  theory  for  subcritical  Reynolds  numbers. 

Region  III  was  highly  turbulent  and  three-dimensional  and  was  very 
difficult  to  interpret  with  certainty.   It  did  appear,  however,  that  a 
separated  region  with  its  characteristic  separation  vortex  was  present. 
No  "trapped"  vortices  were  observed,  but  the  flow  appeared  to  be  somewhat 
stagnant  in  the  region  extending  approximately  30  degrees  forward  of  the 
point  of  attachment  of  the  cylinder  to  the  plate. 
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VI .   CONCLUSIONS  AND  keCdfilfel^DigPIONS 

The  theory  for  the  flow  around  a  right  circular  cylinder  attached  to 
a  flat  plate  in  a  crossflow  of  air  which  was  proposed  in  this  investigation 
appears  to  accurately  describe  the  actual  conditions  existing  at  Reynolds 
numbers  based  on  cylinder  diameter  from  approximately  53,000  to  153,000. 
The  heat  transfer  data  obtained  correlates  well  with  the  pressure  data 
obtained,  and  both  are  in  agreement  with  the  proposed  theory  in  a  qualita- 
tive sense.   Quantitatively,  both  the  heat  transfer  data  and  the  pressure 
data  for  Reynolds  numbers  in  the  range  of  53,000  to  153,000  correlates  in 
the  region  between  0  and  180  degrees  with  the  data  obtained  by  Meyer  [3] 
for  flow  around  a  "free"  right  circular  cylinder  in  a  crossflow  of  air  at 
Reynolds  numbers  in  the  range  300,000  to  500,000. 

The  flow  visualization  exhibited  consistency  with  the  proposed  theory 
for  subcritical  flow,  but  was  inconsistent  in  some  regions  with  the  experi- 
mental results.   This  was  probably  due  to  the  fact  that  the  flow  visualiza- 
tion experiment  was  conducted  at  a  much  lower  Reynolds  number  than  were  the 
heat  transfer  and  pressure  experiments.   It  is  recommended  that  an  attempt 
be  made  to  refine  the  existing  flow  visualization  apparatus  such  that  higher 
Reynolds  number  flows  can  be  attained.   This  would  give  a  clearer  picture 
of  the  flow  at  Reynolds  nuiTibors  in  the  range  of  50,000  to  150,000,  and 
correlation  with  the  other  data  could  be  made  more  meaningful. 

The  liquid  crystal  technique  developed  by  Field  [2]  again  proved  to 
be  an  excellent  method  for  quickly  and  easily  obtaining  heat  trannfer  data 
in  forced  convection  environments.   Since  the  problem  addressed  in  this 
work  had  not  been  investigated  previously,  no  measure'  of  comparison  id 
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available  for  the  values  of  the  local  Nusselt  number  obtained.   But,  the 
qualitative  information  displayed  J?i'~^*^^v.  J. ■'..cp'i.'^.  Ajrvstals  was  verified 
both  by  the  pressure  data  and  the  flow  visualizarxon.   Continued  experiments 
involving  liquid  crystals  are  highly  encouraged. 

The  results  obtained  in  this  investigation  suggest  several  areas  in 
which  follow-up  type  investigations  might  be  conducted.   First,  the  present 
investigation  should  be  extended  over  a  greater  Reynolds  number  range  with 
primary  attention  being  focused  on  Reynolds  numbers  below  50,000.   This  is 
considered  necessary  in  order  to  ascertain  at  what  Reynolds  number  the  flow 
transitions  from  siibcritical  to  critical  for  this  particular  flow  configura- 
tion.  It  is  felt  that  such  a  point  does  exist  because  the  results  of  the 
flow  visualization  conducted  in  this  investigation  suggest  subcritical 
behavior  at  a  Reynolds  number  of  10,000. 

Secondly,  an  effort  should  be  made  to  determine  at  what  distance  away 
from  the  flat  plate  the  cylinder  can  be  moved  and  not  feel  the  effects  of 
the  plate,  i.e.  at  what  distance  away  from  the  plate  will  the  cylinder 
behave  as  a  "free"  cylinder. 

Finally,  a  study  could  be  undertaken  with  the  same  configuration,  but 
with  attention  being  focused  on  the  flat  plate  rather  than  the  cylinder. 
'  This  study  could  be  made  in  the  same  light  as  those  conducted  on  forward- 
facing  and  rear-facing  steps. 


/ 
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HEAT  TRANSFER  DATA  itND  UAin' "KtijruCr ION 

The  values  of  the  local  Nusselt  number  at  various  angular  locations 
were  determined  using  the  relationship  between  the  convective  heat  trans- 
fer coefficient,  h  ,  and  the  surface  heat  flux  produced  by  Joulean  heating, 

c 

V^/RA  ,   This  relationship  was  developed  by  Meyer  [3]  by  performing  a 

s 

simple  energy  balance  on  an  elemental  control  volume  of  the  Temsheet. 
The  result  is 

V^/RA  +  (kt/r  ^)d^T/de 

s _o , 

h  = h 

""  (T-T  ) 


where 

h  =  convective  heat  transfer  coefficient 
c 

V  =  voltage  applied  to  the  test  section 

R  =  electrical  resistance  of  the  test  section 

A  =  surface  area  of  the  test  section 
s 

k  =  thermal  conductivity  of  the  test  section  material 

t  =  thickness  of  the  test  section 

r  =  radius  of  the  cylinder 
0 

T  =  temperature  at  angular  location  9 
T  =  temperature  of  the  air  freestream 

9   =  circumferential  location  on  the  cylinder 

2    2 

h  =  radiation  heat  transfer  coefficient  =  oF        (T+T  ) (T  +T   ) 
r  1-2    '^      °° 

v.'here 

-0     „   ^  2  o  4 
o  ^  Stefan-Boltzman  constant  -  .17] 4  x  10   HTU/hr-ft  -  R 

F    =  radiation  exchange  factor  between  cylinder  and  surroundings 
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Once  h  has  been  determined,  the  Nusselt  nimiber  can  be  obtained  from 
c 

the  equation  '  ^  .'-  ,^  .*.-—- 

Nu  =  2h  r  /k  . 
c  o  air 

The  Reynolds  number  is 

Re  =  U  D/v 

00 

where 

U  =  velocity  of  the  air  freestream 

OO 

D  =  diameter  of  the  cylinder 

V  =  kinematic  viscosity  of  air 

Field  [2]  determined  that  the  product  kt  for  Temsheet  was  on  the  order 
of  0.002  BTU/hr  -   F;  hence,  it  was  considered  negligibly  small  for  purposes 
of  this  investigation. 

A  sample  calculation  of  the  Nusselt  number  is  provided  below  using  the 
data  obtained  at  a  Reynolds  number  of  153,520  and  an  angular  location  of 
330  degrees. 

o    =  .1714  X  10~^  BTU/hr- ft^-°R^ 
T    =  111.0°F  =  571°R 

T     =  73.1°r  =  533. 1°R 

00 

U    •=  76  ft/sec 

OO 

V  =  37.4  volts 

R    =13.1  ohms 

D/v   =  2020  sec/ft 

2r  /k  .   =22.0  hr-°F-ft  /BTU 
o   air 

A    =  0.528  ft^ 

s 


Re    =-•  U^D/v  =  (76  ft/r,ec.)  {?/)?.()   s(-r~/ft)    =  153,520 

.       u   .  ,         V"^  (37..l)^  volts" 

n    =  h  + 


^     ^    RA  (T-T  )     (13.1)  (.528  f  t'^ )  (  37  .  9°F)  ( .  293  watts/BTU/hr ) 
s    <» 
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h  =18.2  BTU/hr-ft^-°F 

h^  =  (.9)  (.1714  X  10~®BTU/hr-ft^-'^R''libTi%-53J-l)  (571^+533.1^)  °R^ 

h  =1.04  BTU/hr-ft^-°F 

r 

h  =17.17  BTU/hr-ft^-°F 
c 

Nu  =  (17.17BTU/hr-ft^-°F) (22.0  hr-ft^-°F/BTU) 
Nu  =  377.7 

The  data  sheets  for  all  nine  heat  transfer  runs  follow. 
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DATA 


Ave.  Air  Temp.  =  61.70P;     D/- =20'/'/ sec 

^  TT 

Air  Speed=.36  cm.H20 

=  aS./  ft/sec 

Re=    SI,  304 

Film  Temp.  =  '7(5.'/ OF;       D/k=}a.l  hr-Op-ft^BTU 


,Date:   10  Stpi  7V 
Time:   /?30 
Resistance=  )3-3j\. 
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DATA 
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DATA 


Ave.  Air  Temp.=  i^Cop;  ■D/f=2.0:\s^ 
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DATA 
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DATA 
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DATA 


Ave.   Air  Teinp.=7'!^5op;     j)A,=^o^Osec 
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APPENDIX  B 
PRESSURE  DATA  AND  DATA  REDUCTION 

The  values  of  the  pressure  coefficient  were  determined  from  the 
formula 

C  =   ® 


P    Js  PU  2 

"      CO 

where 

C  =  pressure  coefficient 
P 

P  =  static  pressure  on  the  cylinder  surface  at  angular  location  9 
6 

P  =  static  pressure  in  the  wind  tunnel  test  section 
s 

p   =  density  of  air 

U  =  air  free  stream  velocity 

2 

Meyer  [3]  had  determined  previously  that  the  quantity   ^j  pU   could 

be  represented  by  the  difference  in  static  pressure  between  the  wind 

tunnel  plenum  chamber  and  the  wind  tunnel  test  section.   Therefore, 

1     2 

^  PU   =  P   -  P 

CO       pc      s 

where 

P   =  static  pressure  in  the  wind  tunnel  plenum  chamber 
pc 

A  sample  calculation  of  the  pressure  coefficient  is  provided  below 

using  data  obtained  at  a  Reynolds  number  of  15G,lfiO  and  an  angular  location 

of  90  degrees 

P    =  12.05  in.  HO 

P    =  4.02  in.  tlO 
pc  2 

P    =  6.58  in.  H^O 
s  2 
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C  =  -<i2.r/^  <-;..6  . 5°  V  /  M .  02  -  6.58) 
P 


C  =  -2.14 
P 


The  data  for  the  two  pressure  runs  follow. 
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PRESSURE  DATA 


Air  Speed=;.5&  cin.H20 
=  51    ft/sec 
Air  Temp.=  6:.r  op 
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PRESSURE  DATA 


Air  Speed=3^  cm.H20 
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Appendix  (-■ 
uncertainty  analysis 


uncertainty  calculations  were  performed  to  ascertain  the  degree  of 
uncertainty  in  the  final  results.   The  method  of  Kline  and  McClintock 
described  in  reference  19  was  used. 

Uncertainties  existed  in  the  following  measured  variables:  voltage, 
resistance,  cylinder  geometry,  surface  temperature,  Temsheet  emissivity, 
angular  location,  manometer  pressure  drop,  all  pressures  in  the  pressure 
experiment,  and  properties  of  air.   The  uncertainties  associated  with  the 
latter  were  considered  negligible. 

TOTAL  HEAT  TRANSFER  COEFFICIENT,  h 

The  total  heat  transfer  coefficient  was  calculated  from  the  equation 

V  /R 

h  = 


A  (T-T  ) 
s     °° 


Therefore,  the  uncertainty  in  the  heat  transfer  coefficient  is 

expressed  as 

I         2       2        2        2 

For  the  experimental  run  at  a  Reynolds  number  of  153,520  and  an  angular 
location  of  330  degrees, 

V   =  37.4  ±  0.1  volts  (20  to  1) 

R  =13.1:^  0.2  ohms  •         (20  to  1) 


o 


AT  =  37.9  ±  0.5   F  (20  to  1) 

A  =  .528  >  0.010  ft^  (20  to  1) 

s 
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r—  =  .027 
H 

h  =  18.2  ±  .5  BTU/hr-ft^-^F       (20  to  1) 
This  result  does  not  take  into  account  the  uncertainty  in  angular  loca- 
tion in  this  region.   This  uncertainty  is  estimated  to  be  ±  3  degrees 
(20  to  1). 

The  uncertainty  in  the  total  heat  transfer  coefficient  remained  rela- 
tively constant  around  the  entire  circumference  of  the  cylinder.   However, 
the  uncertainty  in  angular  location  varied  from  a  maximum  of  ±  5  degrees 
(20  to  1)  in  the  stagnation  region  and  the  wake  region  to  a  minimum  of  ±  1 
degree  (20  to  1)  in  the  regions  near  the  laminar  separation  point  and  the 
reattachment  point. 

RADIATIVE  HEAT  TRANSFER  COEFFICIENT,  h 
r 

The  radiative  heat  transfer  coefficient  was  determined  from  the 

equation 

h  =  oF,  ^(T  +  T  )  (T  +  T   ) 
r     1-2      00        oo 

Assuming  the  uncertainty  in  the  temperature  measurement  is  the  same 

for  both  T  and  T  ,  and  defining  T  =   %  (T  +  T  ) , 

h   =  4aF   T  ^ 
r      1-2  m 

Since  o  is  a  universal  constant,  it  does  not  enter  into  the  uncertainty. 

Therefore, 

r 


h 
r 


The  radiation  exchange  factor,  F    ,  can  be  shov/n  to  be  equal  to  the 
emissivity  cf  the  Tcmsheet.   Since  this  value  wuk  )iot  known,  ci   value  of 
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0.9  was  chosen  based  upon  the  emissivities  of  comparable  materials.   The 
xincertainty  in  this  value  was  assumed  to  be  ±  .05  (20  to  1).   Hence,  at  a 
Reynolds  number  of  153,520  and  an  €ar\C,ui:i:ii-I:v::::a''"  '<  •  of  330  degrees 

\ 

-— ^  =  0.06 
n 

r 

r  =  1.04  ±  0.06  BTU/hr-ft  -°F       (20  to  1) 


NUSSELT  NUMBER,  Nu 

The  Nusselt  number  was  calculated  from  the  equation 

Nu  =  h  D/k  =  (h-h  )D/k 
c         r 

Therefore 


2        2       2         2 

Nu 


Nu 


I)  ■  ft)  •  G)  •  G) 


Since  the  uncertainties  in  the  properties  of  air  were  considered 
neglible  and  the  uncertainty  in  the  diameter  measurement  was  negliglibly 
small. 


Nu 


Nu 


■W^} 


CjJ 

-^     =   .029 
Nu 

Nu  =373  ±  11        (20  to  1) 

REYNOLDS  NUI-SER,  Re 

The  Reynolds  number  was  calculated  from  the  equation 

Re  =  U  D/v 
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Therefore , 


2 
Re 


Re 


y  -0  -Q 


The  quantity  U   is  directly  proportional  to  the  square  root  of  the 
wind  tunnel  micromanometer  reading  AH;  hence 


u  0). 

U  AH 


0.5 


U  AH 

00 

Again,  since  the  uncertainty  in  the  kinematic  viscosity  and  in  the 
cylinder  are  negligible, 

0)  (0. 

_Re  ^   0  5  _AH 
Re     "-^   AH 

At  a  Reynolds  number  of  153,520  (AH  =  3.58  cm.  HO)  the  uncertainty 

in  the  micromanometer  reading  was  estimated  to  be  ±  0.03  cm.  HO.   Hence, 

— —  =   0.5(0.03/3.58) 
Re 

to 

-^  =   0.004 
Re 

Therefore, 

Re  =  153,520  r  G14        (20  to  1) 


PRESSURE  COEFFICIENT,  C 

. . . p 

The  pressure  coefficient  was  determined  from  the  equation 

P„  -  P      AP^ 

c  =  '   "     e 


p    P   -  P       AP 

pc    s       pc 
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Therefore, 


"c 


c 
p 


^   pc/ 


Since  both  the  pressure  differences  were  measured  on  the  same  man- 
ometer bank,  they  shared  the  same  uncertainty  which  was  considered  to  be 
±  0.05  in.  HO  (20  to  1).   Therefore,  at  a  Reynolds  number  of  156,180  and 
an  angular  location  of  90  degrees , 


•^c 

— ^     =   0.022 
C 
P 

C      =   -2.14    ±   0.05 
P 

(20   to  1) 
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